Exact renormalization of the random transverse-field Ising spin chain 
in the strongly ordered and strongly disordered Griffiths phases 
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The real-space renormalization group (RG) treatment of random transverse-field Ising spin chains 
by Fisher (Phys. Rev. 551, 6411 (1995)) has been extended into the strongly ordered and strongly 
disordered Griffiths phases and asymptotically exact results are obtained. In the non-critical region 
the asymmetry of the renormalization of the couplings and the transverse fields is related to a 
non-linear quantum control parameter, A, which is a natural measure of the distance from the 
quantum critical point. A, which is found to stay invariant along the RG trajectories and has 
been expressed by the initial disorder distributions, stands in the singularity exponents of different 
physical quantities (magnetization, susceptibility, specific heat, etc), which are exactly calculated. 
In this way we have observed a weak-universality scenario: the Griffiths-McCoy singularities does 
not depend on the form of the disorder, provided the non-linear quantum control parameter has the 
same value. The exact scaling function of the magnetization with a small applied magnetic field is 
calculated and the critical point magnetization singularity is determined in a simple, direct way. 
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I. INTRODUCTION 

In statistical physics and in the theory of interacting 
many body systems exact solutions are of great impor- 
tance, especially in the vicinity of singular points, such 
as at phase transitions. They provide physical insight 
about the cooperative processes and their results could be 
used as testing ground for different approximations and 
numerical methods. Exact solutions about models with 
short range interactions and in the presence of quenched 
randomness are scarce which has greatly hampered our 
understanding about collective phenomena in disordered 
systems. At present remarkable examples about exactly 
soluble problems in the above class are the criticaLbehav- 
ior of low-dimensional random quantum systemsu. Here 
the interplay of quenched disorder, quantum fluctuations 
and correlations can be systematically studied within a 
real-space renormalization group (RG) scheme, which is 
expected to lead to asymptotically exact results, at least 
for strong enough disorder. The prototype of such type of 
random quantum systems is the random transverse-field 
Ising spin chain (RTIC) for which perhaps the most de- 
tailed analytical and numerical information is available, 
as far as the random quantum critical behavior is con- 
cerned. The RTIC is defined by the Hamiltonian: 



H 



a, cr 



i+1 



(i.i) 



in terms of the af' z Pauli matrices at site, I, whereas 
the transverse fields, hi > 0, and the couplings, J; > 0, 
are independent random variables taken from the (ini- 
tial) distributions, P m {h)dh and R m (J)dJ, respectively. 
The (linear) quantum control parameter of the model is 
defined asB 



[ln/i] av -[lnJ] in- 
var [In h] + var [In J] 



(1.2) 



where var(x) is the variance of x and here and in the fol- 
lowing we use [. . .] av to denote averaging over quenched 
disorder. The quantum critical point at S = separates 
the ferromagnetic (6 < 0) and the paramagnetic (5 > 0) 
phases. 

Some previously known exact results about the surface 
magnetizatioroQ and about tynical correlations!^ have 
been greatly extended by FisherO using the strong disor- 
der RGr*nethod, originally introduced by Ma, Dasgupta 
and HuQ. Fisher has shown that at the critical point 
the distribution function of the couplings and that of the 
transverse fields broaden without limits as the energy 
scale, f2, as defined by the strongest bond or transverse 
field, is lowered. Therefore, as the fixed point of the 
transformation with = is approached the disorder 
becomes stronger and stronger, so that in this, so called, 
infinite randomness fixed point (IRFP) the ratio of typi- 
cal couplings and transverse fields at neighboring sites is 
either zero or tends to infinity. As a consequence in the 
IRFP the RG transformation becomes asymptotically ex- 
act and the fixed-point RG equations for the RTIC can be 
solved, in large extent, analytically. From the RG treat- 
ment and from other analytical resultsu we have a clear 
physical picture about the origin of the critical behavior 
of the RTIC, which is most probably also relevant for 
other low-dimensional random quantum critical systems, 
having an IRFP. The critical properties of these sys- 
tems are determined by the so called rare events, which 
are realizations occouring with a vanishing probability, 
but which dominate the average properties. In contrary, 
the typical realizations, which appear with probability 
one have a vanishing contribution to the average critical 
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quantities. The mathematical origin of the solubility of 
these models is connected with the above observation, 
since it is enough only to deal with the rare events and 
an overwhelming part of the realizations are irrelevant in 
respect of the average critical properties. 

In the following we list the existing results about the 
singular behavior of the RTIC. At T = and with- 
out an external field, H — 0, the average magnetiza- 
tion of the surface spin, m s , asymptotically vanishes as 
m s ~ (— 5)P S , with an exponent! 



0s = 1 



(1.3) 



whereas the same behavior fox the bulk magnetization, 
to, involves the bulk exponenta, (3: 



/3 = 2-$, $ = 



1 + V5 



(1.4) 



Average correlations, G(r) = [(erf of^ r )] av , outside the 
critical point decay exponentially and the correlation 
length, £, asymptotically diverges in the vicinity of the 
critical point asQ 



(1.5) 



On the other hand the typical correlation length, ^typ, 
defined through the relation [In ((erf crf +r ))] av ~ ~r/£ tyP 
involves another exponcntEl, v tyv = 1. At the critical 
point average correlations decay as a power, G(r) ~ r~ 2x , 
and the scaling dimension of the bulk magnetization, x, 
satisfies the scaling relation x = (3/v. Similarly, for end- 
to-end critical correlations the corresponding scaling di- 
mension, x s , is expressed by x s — f3 s /v. 

In a quantum system statical and dynamical correla- 
tions are inherently related. In the RTIC at the critical 
point dynamical scaling is strongly anisotropic, the rel- 
evant (imaginary) time scale the relaxation time, t r , is 
related by the length-scale, £, as 



In t r 



1/2 



(1.6) 



In the RG study Fisher has also considered the weakly 
ordered and weakly disordered Griffiths phasesa, which 
are situated in the vicinity of the critical point. He found 
an anisotropic scaling relation, t r ~ ^ 1 / 2 \ s \ j which can 
be interpreted as a 6 dependent dynamical exponent, z, 
which in leading order is given by 



1 

W\' 



\s\ « i 



(1.7) 



In the presence of a small external field, H <C 1, and 
in the vicinity of the critical point, \8\ <C 1, but with a 
finite combination of, 7 = 5\n(H) = 0(1), Fisher has 
obtained the exact scaling function of the magnetization 
of the form: 

m(6,H) » m MH /H)f~ 2 m[6 HH /H)} , (1.8) 



where too and Hq are dimensional constants, and 771(7), 
is given in Ref. |^. 

In the present paper we extend the RG treatment by 
Fisher into the entire Griffiths region, which for some 
type of initial disorder distribution could cover the entire 
off-critical region, < \6\ < 00. Our analytical results in 
the Griffiths phases are asymptotically exact in the same 
sense as argued by Fisher at the critical point. Here we 
summarize our main findings. A short account of our re- 
sults, especially about renormalization of couplings and 
transverse fields, has been announced in Ref. |l^. 

The asymmetry in the renormalization of the couplings 
and the transverse fields is related to a non-linear quan- 
tum control parameter, A, which is defined as the root 
of the following equation: 



= 1 



(1.9) 



We have shown that A stays invariant along the RG tra- 
jectories (i.e. as the energy scale, f2, is lowered) and so 
can be expressed by initial disorder distributions. Note 
that for very different disorder distributions one might 
have the same A, provided the distributions have the 
same form of asymmetry. The dynamical exponent, z, 
is simply expressed by the non-linear quantum control 
parameter as: 



1 
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(1.10) 



which is an exact relation in the entire Griffiths region. 
We have calculated the singular behavior of different 
physical quantities (magnetization, susceptibility, spe- 
cific heat, etc) and the singularity exponents are all ex- 
pressed by the non- linear quantum control parameter, A. 
In this way we have demonstrated a weak-universality 
scenario: details of the disorder distributions are ir- 
relevant for the Griffiths-McCoy singularities, provided 
the non-linear quantum control parameter has the same 
value. 

In the vicinity of the random quantum critical point 
the non-linear and the linear quantum control parame- 
ters are asymptotically identical, 



A = S + 0(S 2 ) 



(1.11) 



consequently in the linear A limit from our formulae we 
can recover Fisher's results^ about the weakly disordered 
and weakly ordered Griffiths phases. From our results one 
can also obtain the scaling function of the magnetization 
as a function of a small applied magnetic field, H, which 
in the paramagnetic phase is given in the functional form: 



m(A,H)=m D ( — 



2A 



A 

Ad ' 



H 



2A 



(1.12) 



where tod, Hd and Ad are non- universal dimensional 
constants. Here Ad is a non- universal parameter, which 
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depends on the details of the distribution of th e diso rder. 
At the critical point, A = 0, o ur result in Eg. ( 1.12 ) goes 



over Fisher's result in Eq.([L8|). 

The structure of the paper is the following. The RG 
equations and their solution for the fixed point distribu- 
tion of the couplings and the transverse fields are pre- 
sented in Sec. 2. Renormalization of lengths and mag- 
netic moments are given in Sec. 3. The scaling behavior 
of different thermodynamic quantities in the presence of 
a finite external magnetic field or at small, but non-zero 
temperature is calculated in Sec. 4. We conclude our 
paper with a discussion about possible extension of our 
results to other problems in Sec. 5. Some detailed cal- 
culations about the distribution function of lengths are 
given in the Appendix. 



II. RENORMALIZATION OF COUPLINGS AND 
TRANSVERSE FIELDS 

Here we co nsid er the RTIC as given by the Hamil- 
tonian in Eq.(l.l) without external magnetic field, i.e. 
H = 0. To a spin at lattice site, I, we assign a magnetic 



moment, /i;, and a length, If, while the transverse field 
acting on this spin is denoted by hi. Similarly, to the l-ih 
bond, connecting lattice sites I and I + 1, we assign the 
length, l\, and the associated coupling is denoted by Ji. 
In the initial situation If = If = 1/2, fii = 1 and the cou- 
plings and fields are taken from the initial distributions, 
Pin(h)dh and i?; n (J)dJ, respectively. 

During renormalization the strongest term in the 
Hamiltonian, coupling or transverse-field, of strength 
f2 is successively decimated out and the neighboring 
transverse-fields or couplings are replaced by weaker 
ones, which are generated by a second order perturba- 
tion calculation. If the strongest term is a coupling, say 
Q = Ji, then the two spins connected by J; flip coher- 
ently in a longitudinal field, thus they act as an effective, 
composite spin having the renormalized parameters: 

fc=^±i, V' = l! + il + l!+ 1 , a = w + w+i ■ (2-1) 

On the other hand, if the strongest term in the Hamil- 
tonian is a transverse-field, say Vt = hi, then the state 
of this spin in a small longitudinal field is practically 
unchanged, thus its contribution to the susceptibility is 
negligible. Consequently in magnetic point of view this 
spin can be decimated out, and the renormalized param- 
eters of the new effective bond connecting sites I — 1 and 
I + 1 are given by: 



J 



<h-i<h 
~hT' 



I s = lit + It + li 



(2.2) 



Note that the decimation equations for h and J are re- 
lated through duality. 

During renormalization the energy scale is reduced 
and the joint distribution functions, such as for the 



spins, P(h, I s , fx; Q,)dhdl s dfi, and that for the bonds, 
Ri(J,l h ;Sl,)dJdl b , are also fi dependent. Generally we 
deal with the following reduced distribution functions: 

P (h,Q) = J J P{h,l s ,(i;Sl)dl s dn 

Pi{Ki s ,n) = J p(h,i',msi)dn 
p»{h,n,n) = J p(h,i s ,fi;n)di s 

Ro(J,n) = J Ri{J,l b ; n)dl b , (2.3) 

all of which are normalized. In this Section we consider 
the distribution of transverse fields and couplings, so that 
we work with Po(h,Q) and Ro(J,Q,), whereas the other 
joint distributions, which are connected to the size of av- 
erage lengths and average moments will be considered in 
the following Section. 

We start to calculate the variation of the distribution 
function of transverse fields, dPoQi, fi), when the energy 
scale is lowered by Q — > O — df2, which amounts to elim- 
inate a fraction of df2[Po(fi,f2) + Ro(£l, f2)] spins. Here 
one should take into account the fact that as a strong 
bond is decimated out two original fields are also elim- 
inated and one new is created, the strength of which is 
given in Eq. fl2.ip . Since the new distribution function 
should also be normalized we arrive to the equation: 



d P 

= P (h, Q) [Ro(Q, n) - P (O, fi)] 

iJo(n,fi) / dh'p {h',n)p (^n,n)^j . 



(2.4) 



One can similarly derive the evaluation equation of the 
coupling distribution: 



dn 



Ro(j,Q) [p (n,n)~R Q (n,n)} 



-p (n,si)J dJ'R (j',n)Ro(^n 7 n)y , (2.5) 



which follows simply from Eq.(2.4) by duality, which 
amounts to interchange h «-> J and Pp < -> Rq. Th e two 
integro-differential equations in Eqs.(2.4) and (2.5) have 
to be supplemented by the initial conditions, represented 
by the distributions Pi n (h) and i?i„(J). 



A. Fixed-point solution 



A special solution to the problem in Eqs.(2.4) and (2.5) 
is given by the functions: 



p (h,n) 



n \h 



(2.6) 
(2.7) 
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thus they depend only on the values of the distributions 
at their edges, at P (Q,£l) = Po/Q and at R (Q,£l) = 
ro/fi. At the end of the Section we present arguments, 
that this special solution represents the true solution of 
the problem at the fixed point, i.e. as fl — > 0. Later 
we also show how the parameters of the special solution 
can be related with the initial distributions, Pi n (h) and 
i?i„(J). 

Putting Eqs.fp^) and (O) into Eq.(E3) we obtain: 



p Q r Q - SI 



dpo 
dfl 



ln- 



fi 



1 

Po 



Q 



(2- 



which leads to the ordinary differential equation: 

-Poro , (2.9) 



dp 

dr 



in terms of the log-energy variable, T — — In fi. Similarly, 
from Eq.(2.5) we obtain for the edge-parameter, rp: 



dr 

dr 



-r Q po . 



(2-10) 



Subtracting Eq.(2.9) from Eq.(2.1C) we obtain that po 
and ro differ from each other by a constant, 2A, 



Po - r = 2A . 
Thus in terms of the variable 

Vo = Po - A = r + A 
we obtain one differential equation: 

A 2 . 



dyo 

dr 



(2.H) 



(2-12) 



(2-13) 



Here we note that A is related to the asymmetry in the 
renormalization of couplings and transverse fields and its 
value, which can be expressed by the initial distributions, 
will be determined later. At the critical point, where the 
distributions of the transverse fields and that of the cou- 
plings evolve to the same limiting function as fi — ► 0, we 
have A = 0. In the paramagnetic phase, where according 
to Eq.(1.2), transverse fields in average are stronger than 
the average couplings, A > 0, whereas in the ferromag- 
netic phase we have the opposite situation, A < 0. 



At the critical point, A = 0, the solution to Eq.(2.13) 
is given by: 



1 



1 



Vo=Po = r a 



T-T Q ln(fi /fi) 



5 = A = 0, 



(2-14) 



where To = — In Qq is a reference (log)energy scale. 
It is instructive to consider the distribution of the re- 
duced log-coupling variable rj = — (lnfi — In h ) / In fi = 
— (lap, — In J)/ lnfi, which is given from Eqs. fl2.6[ ) and 
(PH ) as 



p(rj)dr] = exp(— r])dr] 



(2-15) 



This is. just the critical point solution of the RTIC by 
Fisheitl 

The solution to Eq. fl2.13 ) in the off- critical region, 
A 7^ 0, is given by: 



Vo 



Ay + A 2 th [A(r - r )] 
A + y th[A(r-r )] 

|AI j 1 + 2^4 (W)) 2A 



(2.16) 



where the solution goes through the point yo = y at the 
refe rence (log) energy cut-off, IV The second equation in 
Eq.( 2.16] ) is the approximate form of the solution close 
to the line of fixed points, where in terms of the original 
energy-scale variable fi/fi <C 1. We note that yo/A is a 
unique function of two dimensionless variables 



yo 
A 



A 



fi 



2:S 



(2-17) 



where Ap = yo and Qd = fio- 

In the following we relate the asymmetry parameter, 
A, with the properties of the initial distributions, P m (h) 
and R in (J). For this purpose we calculate first the deriva- 
tive: 



d 



dRo(J,n) 
dfi 



J^dJ 



[J»} av = R (Sl,n)W+ f 
Jo 

= R (n, + (p (fi, fi) - R Q (n, fi)) [J"] av 
-p (fi,fi) J djr J dj'R(—n,n\ — , (2. 



18) 



where in the second equation we have used the RG equa- 
tion in Eq.( |2.5| ). In the last term we change the order of 
the integration and thus obtain: 



jfi 



dj' / djj^R — fi, fi — = n-v [j^] 



] 2 

J av 



(2-19) 



One can evaluate in a similar way (d/dil)[/i _A1 ] av and 
then obtain for the average value of the following deriva- 
tive: 



_d_ 

dfi 



h 



= 1- 



h 



,} . (2.20) 



Notice that this quantity is vanishing for a parameter 
(i = p,, provided [(J/h)^] = 1, thus fi defined in this 
way stays invariant along the RG trajectory. This rela- 
tion is valid even if the starting RG steps are approxi- 
mative. At the fix ed point, fi — > 0, with the solution in 
Eqs.([D]) and (p^l) we have evaluated the average: 



4 



= 1 



(2.21) 



Consequently the asymmetry parameter, A, is an invari- 
ant quantity of the RG transformation, which is then 
dete rmi ned by the initial distributions, as announced in 
Eq.(1.9). In the following we call A as the non-linear 
quantum control parameter of the RTIC. 



Next we show that the RG equations in Eqs.(2.1) and 



(2.2) become asymptotically exact as the line of fixed 
points is approached, i.e., as O/Oq ~^ 0. Let us here 
consider the disordered Griffiths phase, A > 0, the rea- 
soning for A > follows from duality. Here the ratio of 
decimated bonds, Anj, and decimated transverse fields, 
An/,, goes to zero as Anj/Arih — i?o(ft, ft)/-Po(ft, ft) = 
fo/po ~ ft 2A , thus close to the fixed point almost ex- 
clusively transverse fields are decimated out. Then the 
probability, Pr(a), that the value of a coupling, J, being 
neighbor to a decimated transverse field is ft > J > aft 
with < a < 1 is given by 

Pr{a) ~ / Rq(J, ft)dj = 1 - a r ° « r ln(l/a) , 
Jan 

(2.22) 

whic h goe s to ze ro d uring iteration, since according to 
Eqs.(^16|) and ( [H^ r = Ro(Q,Q)Q -> 0. Conse- 
quently the RG transformation becomes asymptotically 
exact and the singularities, which are characterized by 
the paramet er A as calculated by the original distribu- 
tions in Eq.(2.21), are also exact. 



B. Relation between energy- and length-scale 

Next we are going to study the actual relation be- 
tween the asymmetry or non-linear quantum control pa- 
rameter, A, and the Griffiths-McCoy singularities of the 
RTIC. For this we investigate the relation between the 
energy scale, ft, and the length scale, Lq, by studying 
the fraction of non-decimated spins, na- When the en- 
ergy scale is decreased by an amount of dft a fraction of 
spins, dnn = nn[Po{Q, ft) + Ro{Q, ft)], is decimated out, 
so that we obtain the differential equation: 



dft 



= n fi [P Q (ft,ft)+.R Q (ft,ft)] , 



(2.23) 



what can be rewritten as 
din nn 



dm ft 



(r (n)+p (n)) = -2 W ,(n). (2.24) 



Using the solution to yo(fl) in Eq.(2.16) one can integrate 
Eq.Q2.24p with the result: 



ch 



A In 



fto 

ft 



A 



A In 



ftp_ 
ft 



(2.25) 



To obtain the scaling form at the critical point we take in 
Eq. fl2.25D the limits, A -> and r = - In ft -> oo, with 
however Axf->0 and obtain: 



no 



1 + 2/ In 



fto 
ft 



In 



fto 

ft 



A = 0, (2.26) 



which could be also directly calculat ed fr om Eq.(2.25) 
wit h the critical point solution in Eq.( 2.14 ). Thus, from 
Eq.(2.26) we get for the typical distance between remain- 
ing spins, L n , as: 



L a 



In?? 
ft 



A = 0, 



(2.27) 



which js just the relation in Eq.( |l.6| ) as found earlier by 
FisheiO. 

I n the Griffiths phases, |A| > 0, one obtains in 
Eq.( |2.25D , n n ~ ft 2 l A l, in the limit ft -> 0. Consequently 
the relation between typical distance between remaining 
spins, Ln ~ l/nn, and the energy scale is given by: 



WA + z/o)' 



fto 
ft 



2|A 



ft~ 2|A| . (2.28) 



Thus A is simply related to the dynamical exponent, z, 

1 



2IAI 



(2.29) 



as announced in Eq.( [1.10 ). 

To obtain a relation betwee n the non-linear quantum 
control parameter, A. in Eq.( 2.21 ), and the linear con- 
trol parameter, <5, as defined in Eq.Jl.2D we perform a 
Taylor-expansion for: 

[J^] av = 1 + 2A [In J] av + 2A 2 [In J) 2 ] av + 0(A 3 ) , 

(2.30) 



and similarly for [h 2A ] 
we obtain, that A(5) 



Putting these into Eq. fl2.2l[ ) 
S + 0(5 2 ), as announced in 
Eq.fll.llD. 

Closing th is s ection we argue that the special solu- 
tion in Eqs.( ^.6| ) and ( |2.7| ) is the true fixed-point solu- 
tion as ft — > 0. n First, we refer to Fisher's results at 
the critical pointO, which justifies that any non-singular 
initial d istr ibution is attracted by the special solution 
in Eqs.(2.6) and (2.7). Second, we consider a finite 
energy- scal e, ft > , when the non-asymptotic solution 
to E cis.(|2.4 |) an d (|2.5| ) are given by the special solutions in 
Eqs.( |2.6|) and ( [2.7P extended by non-universal functions, 
PqQi, ft) and R' (J, ft), respec tive ly. Inserting these non- 



asymptotic solutions into Eq. (|2.6|) the relation in Eq. ( 2 



will be extended by other terms, containing Pg and R' n 



As ft goes to zero, however, the second factor in Eq.(2 
is diverging, therefore t he c orrections become irrelevant 
and the relation in Eq.(|2.9|) will govern the fixed-point 
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behavior. Our thir d arg ume nt is based on numerical so- 
lutions of Eqs.(2.4) and (2.5), wh ich are evo lving towards 
the special solutions in Eqs.(2.6) and (2.7) for different 
initial distributions. 

Thus we can summarize that in the Griffiths phases 
the RTIC is uniquely characterized by a non-universal 
quantum control parameter, A, whi ch is related to the 
dynamical exponent, z, through Eq.( l.lO ). The possible 
difference between two initial distributions leading to the 
same A is given by the non- universal parameters, fio and 
y , which account for the number of neces sary RG s teps 
until the fixed point distributions in Eqs.(2.6) and (2.7) 
are sufficiently approached. 



III. RENORMALIZATION OF LENGTHS AND 
MAGNETIC MOMENTS 

The scaling behavior of lengths and magnetic mo- 
ments during renormalization can be deduced from the 
joint distribution functions, Pi{h, I, fi), Ri( J, Z, fi) and 
Pfj,(h, /i, fi), as defined in Eqs.(2.3). From here on we 
drop the index s or b to indicate the type of the length. 



A. Scaling of lengths 

In the following we consider the joint distribution, 
Pi(h, I, fi), and write down the relevant evaluation equa- 
tion when energy scale is lowered as fi — > fi — dSl. Gen- 
eralizing the reasoning leading to Eq.(2.4) we obtain: 



<LFj(M,n) 

dfi 



= p,(h,i,n) [Ro(n,n)-p (n,n)] 



n 



1-1-2 



dli-R,(fi, I 2 , fi)A(/n, ii.fi) 



dh x — / d/ 2 

h "i Jo Jo 



xPi ^—n,i-h-i 2 ,nj , (3.i) 

and similarly for the coupling distrib utio n by interchang- 
ing R <-* P and J t <-> hi, as in Eqs.(fO|) and (2.5). The 
second term in the r.h.s. of Eq.(3.I) can be written as a 
convolution in terms of the variable I' = h + I3 = I — 1% 

as 



- f dl'Ri(fl,l-l',fl) [ dh 
Jo Jo 

x f dh£pi{hi t h,n)Pi£si,r -h,si) . 

Jh h i h i 
Consequently taking the Laplace-transform: 



(3.2) 



Pi(h,X,ft) = / e- lx Pi(h,l,fl)dl 
Jo 



e- lx Ri(J,l,Vl)dl , 



(3.3) 



we obtain a simpler relation 
dPi(M,fi) 



dil 



Pi(h,X,Q) i?j (fi, 0,fi)-fi(fi, A, fi) 



^ , h 



fi 



-R l (n,x,vi)J^ dh'Pi(h',x,n)Pi(-n 7 x,rt)- , (3.4) 

and similarly for the coupling distribution: 

dRi(j,x,rt) 



dn 



Rt(J,\,Q) P ; (fi,0,fi) - Ri(Cl, A,fi) 



-Pi((i,\,si)J dJ'R i (j',x,n)R i (ySi,x,n)y . (3.5) 

Note that the different A components are separated, 
which makes possible to solve the equations. For A = 0, 
when Pi(h,0,Sl) = P {h,n) and i? ; (J,0,fi) = i? (J,fi), 
the solutions are given in Eqs.(p| and With this 

guidance we are looking for the solution for general A in 
the form: 



Pi(h,x,n) 



pi(X,ft) /fi 



fi 



j 



l- Pl (\,Q) 



i-n(A,n) 



(3.6) 
(3.7) 



where now pi(0, fi) = 7q(0, fi) = po(Cl) and n(0, fi) = 
pi(0,Q) = ro(fi), whereas for A > pt (X, fi) > 7r ; (A, fi) 
and r;(A, fi) > pi{X, fi). This latter relation follows from 
the fact that the average length of a bond, \ > 0, and 
that of a spin cluster, 7 S > 0, is given by: 



h = lim - 

a^o A 



MA,fi) 



rz(A,fi) 



l s = lim i 

a^o A 



f - 



^(A,fi) 
Pi (A, fi) 
(3.8) 



Inserting the functions in Eqs.(3.6) and (3.7) into 
Eqs.(p.4|) and ( p^ we obtain a set of ordinary differential 
equations: 



dpi 

dr 

dn 

df 



dn l 

dr 

dpi 

dr 



-TTl(pi - TTl + pi) 

-pi( r i - Pi + n) ■ 



(3.9) 



involving the functions pi, ni , ri and pi . 

These equations are solved in the Appendix. Here we 
only consider the scaling behavior of the average lengths 
and for this it is enough to treat the small A expansions 
up to linear order: 

Pi(X,n) =p (fi) + Api(fi), 7T*(A,fi)=p (fi) + A7ri(fi) 
n (A, fi) = r (fi) + An (fi), pi ( A, fi) = r (fi) + X Pl (fi) , 

(3.10) 



Inserting the expressions in Eqs.(3.10) into Eq.(3.9) we 
obtain for the correction terms: 



G 



dpi d7T! 

— = -7Tir - pipo, -7j=T = -TTiro -po(Pi - TTl + Pi) 
dr-i dpi 

— = -7Tir - pip , -^r = -PlPo - M r i - Pi + TTl) • 

(3.11) 



Now noticing that 
dln(pi - 7Ti) 

dr 



dln(ri - pi) 

dr 



(3.12) 



we obtain after integration using Eqs ■ ( J2 . 9| ) and ( 2.10 ), 
that pi — 7Ti = Ap/ro and ri — p\ = A r /po, where 
A p , A r are integration constants. Consequently the aver- 
age lengths from Eq.(3.£) are given by: 



L = 



and 



Pi ~ 
Po 



Z s (fi 



r (fio)po(fio) 



n - pi 

^0 



Z 6 (fi ) 



ro(fio)po(fi ) 
r (fi)p (fi) 



Vo 



^(fio)^§ 



-A 2 ' 
(3.13) 

-A 2 



vi - a 2 • 

(3.14) 



At the line of fixed points, fi — > 0, one can see that 
Is ~ h ~ consequently the prev ious interpretation of 
the dynamical exponent, z, in Eq.(2.2S) is justified also 
with the average lengths-scales. 

Now to calculate the correlation length. £, in the para- 
magnetic phase, A > 0, one should take into account 
that the ratio of (non-decimated transverse-fields) / (non- 
decimated couplings) at an energy-scale, fi, is given by 
Po/fo- Consequently the number of non-decimated spins 
in a cluster is given by ~ T s po/ r o ~ l/ r o ~ A -2 , which 
stays constant as the energy-scale is lowered. This quan- 
tity is actually the measure of the size of the average 
correlated domain in the system, where the couplings be- 
tween the spins, being larger then the transverse fields, 
are decimated out. Therefore in this way we have an 
estimate for the correlation length close to the critical 
point: 



(3.15) 



which is consistent with Fisher's result in Eq.(1.5) 



B. Scaling of magnetization moments 

In this subsection we perform a similar calculation 
about the joint distribution function, P^(h, fi, fi), and 
calculate the average size of a magnetic moment, /l(fi), 
as a function of the energy cut-off. The joint distribution 
function, P^h, p, fi), satisfies the differential equation: 



dP„(fe,/z,n) 

dfi 



= p ll (h,n,si) [i*o(n,ft)-p (n,fi)] 



Ro(n,Si)J dh'^ df J ,'P IM (h', t i',Q)P(^Q,fi-n',il) , 

(3.16) 



what can be derived along the lines of Eqs. (2. 4 ) and ( p.l| ). 
The second term in the r.h.s. of Eq.( |3.16 ) is a convolu- 
tion therefore we introduce the Laplace-transform: 



e-^P^p, fi)dp 



(3.17) 



which satisfies the relation 

dP M (M,n) 



dfi 



= p„(M,n) R (n,n)~p^(n,o,n) 



h 



n 



Ro(Q,n)J dh'p^h', s ,n)p^-n, s ,n)- . (3.18) 



In Eq.( 3.1£ ) the different s components are separated, for 
s = 0, when P„(/i, 0,fi) = Po(/i, fi) the solution is given 
in Eqs. (2. 6). As for the joint distribution of the lengths, 
Pl(h, A, fi), in Eq.(p.(^) we are looking for the solution for 
general s in the form: 



P fi (h,s,Q) = 



7r M (a, fi) ^fi 

fi - 



1— P^(s,f2) 



(3.19) 



Here again p jU (0, 0) = ^(0,0) = po(fi) 5 whereas 
p /i (s, fi) > 7r M (s,fi) for s > 0, since the average cluster 
moment, p > 0, is given by: 



p = lim — 

s^O s 



Pti{s,^) 



(3.20) 



Putting Eq. (|3.19| ) into Eq. (p^ ) we find that the func- 
tions, p M and 7r M , satisfy the differential equations: 



= -itn,r 



dp^ 

dr 

d*V / i \ 

— = -7i>(r -p +Pn) 



(3.21) 



Keeping in mind that the average cluster moment, p, and 
thus the average magnetization, to, defined as: 



Is 



(3.22) 



is re lated to the small s asymptotics of the distribution in 
Eq.( 3.19 ), we perform the expansions up to linear order: 



P AI (s,fi) =p (fi) +SPl(fi), 7T M (s,fi) =p (fi) +S7Ti(fi) . 

(3.23) 

For the correction terms, pi and tt i , we derive differential 
equations in terms of the functions po and tq as: 



7 



dpi 



which leads to: 



(J/0 



7Tl 
PO 

Po 

2^ 



(3.24) 



Pi 



Pi 



(3.25) 



where y$ = yo(Q) is given in Eq.([OE). We note that 
with the solution for p\ we have for the average cluster 
moment: 



l" = 



P i ~ 7Ti 
PO 



4 w °d2/oPi(yo)/(yo - A ) 



yo + a 



(3.26) 



ylt i/ie critical point with A = the solution of Eq.(3.25) 
is given in a simple power-law form: 



Pi = Vo 



0, 



(3.27) 



where r = (V5— 1) /2 is the positive root of the equation: 
t(t + 1) = 1. (The other linearly independent solution 
with t = — (v5 + l)/2 is physically not acceptable, since 
the average cluster moment would be smaller than 1.) 



From Eq.( |3.25| ) we have tv\ = —ry^ 7 and using Eq.( 3.26 ) 
we obtain for the average cluster moment at the critical 
point: 



\x = const y 



Mo 



hi 



1 + V5 
2 

(3.28) 



In this way we have re-derived Fisher's results about the 
scaling behavior of the average cluster moment in a direct 
way. 

In the Griffiths phases with A ^ the differential 
equation in Eq.( 3.25| ) in terms of the variable y = yo/A 
is related to the Legendre differential equation and the 
physically acceptable solution can be expressed by the 
hypergeometric functionO, F(a, b; c; z), as 



Pi = |A| 



-^1-,- + -;- 

2' 2 2' 2 



t\\ I = A! 

y 



7i(y) , 

(3.29) 



where, in the limit A —> we recover the solution at the 



critical point in Eq.( 3.27 ). From Eq.( 3.25 ) we obtain 



IA 



<t>i{y) , 



1 

' 2 

V 

(3.30) 



and putting Eqs. fl3.29| ) and ( |3.3C| ) into Eq.( |3.26| ) we ob- 
tain for the average cluster moment 



fj, = const | A | 



-,-x M 

y + l 



(3.31) 



where f(y) = fi(y) — 4>i{y)- Here one should differenti- 
ate between the paramagnetic (A > 0, y > 0) and the 
ferromagnetic (A > 0, y > 0) phases. In the former case 
the average cluster moment is approaching a finite lim- 
iting value, as Q/Qq — > 0, whereas in the ferromagnetic 
phase, where y — * 1~ in the fixed point, thus /Z is diver- 
gent, as ~fi{£ l) ~ £! -2 l A l. For the average magnetization 
in Eq.(3.22) then one obtains: 



m = mo 



a-y)f(y) 
d-y)f(y) ' 



(3.32) 



where mo is the average magnetization at f2 = Qq 
and y denotes the value of the variable y at the same 
energy-scale. The average magnetization in the param- 
agnetic phase is zero, wher eas in the ferromagnetic phase 
one have to evaluate Eq.( 3.32 ) along the line of semi- 
critical fixed points, O/Oq - ¥ 0. Here taking the limit 
|A| <C 1, i.e. being close to the criticaLpoint we have 
(1 - y)- 1 ~ |A| and f{y) - |A| T so thaM 



in 



constlAI 1- 



const|<5| 1_ 



(3.33) 



From Eq.jOal) one can read the critical exponent of the 
average magnetization as: 

/3=1-t = 2- $, (3.34) 

which corresponds to Fisher's result in Eq.( |1.4| ). 



IV. SCALING OF THERMODYNAMIC AL 
QUANTITIES 

In the previous Sections we have presented the solution 
of the RG equations in the entire Griffiths region for the 
distribution of couplings, transverse fields, lengths and 
magnetization moments. Then, with those distributions, 
average quantities, such as length-scales, magnetization, 
etc. are calculated at zero temperature and in the ab- 
sence of a longitudinal magnetic field. In this Section we 
extend these calculations and determine the scaling form 
of singular thermodynamic quantities as a function of a 
small, but finite temperature, T > 0, or magnetic field, 
H > 0. 

To treat the effect of a small finite temperature in the 
RG scheme one should first notice that the thermal en- 
ergy sets in an energy scale, Ot ~ T, and the RG deci- 
mation should be stopped as is lowered to f^T- At that 



energy scale a fraction of spin clusters, n^ T , in Eq.(2.25) 
is not decimated out and these spins are loosely coupled 
comparing with the temperature, T. Consequently the 
entropy per spin, s, is given as the contribution of non- 
interacting spin clusters: 



s ~ n nT In 2 



(4.1) 



S 



whereas the specific heat can be obtai ned throug h de riva- 
tion: c v = T(ds/dT). From Eqs.(p]) and ( (2~25[ ) we 
obtain for the singular behavior: 



and similarly for the susceptibility: 



s(T) ~ c v (T) ~ T 2|A| 



(4.2) 



which is valid both in the ordered and in the disordered 
Griffiths phases. 

Next, we consider the effect of a small longitudinal 
field, H > 0, at zero temperature. During rcnormaliza- 
tion the local longitudinal field, Hi, at site I is trans- 
formed as 



Hi = Hm , 



(4.3) 



so that the energy-scale related to the longitudinal field 
is given by VIr — H~p(Q). As il is lowered to Qh, i.e. 
when the energy scale satisfies the equation 



(4.4) 



the RG procedure is stopped and the remaining spin clus- 
ters are practically uncoupled. Then the average magne- 
tization and the average susceptibility satisfy the equa- 
tions: 



m(H)=m(a = fl H ), X =^ 



(4.5) 



In the disordered Griffiths phase, where ~P(Qh) has a SIh 
independent limi ting value, we have ~ H, conse- 
quently from Eq.( 3.32| ) the singular behavior is given by 



2A 



A > 



(4.6) 



More generally the scaling form is given in Eq.(1.12), 
where the scaling function can be computed using 
Eqs.(|3~32l), ( pTl6| ) and ( p7| ). Similarly one obtains for 
the scaling of the susceptibility in the disordered Griffiths 
phase: 



-1+2A 



X (T)~T- 1+2A , A>0, 
(4.7) 



where the temperature dependence follows from the scal- 
ing relation, ~ ^t- 



In the ordered Griffiths phase, where ~ ^ 



-2|A| 
H ' 

as given above Eq.( |3.32|) we have n H ~ £f 1 /( 1 + 2 l A l). 
Putting this result into Eq.( |3.32 ) and using the asaimp- 
totic expa nsion for the hypergeometric functions!!!! in 
Eqs.(3.29) and (3.3C) we obtain for the leading field de- 
pendence of the magnetization: 



H \ 1 + 2|A| / JJ 

,»(//! (— J lnf — 



A < 



X(H)~(-) ln(-j, A<0. (4.9) 

Note that in the ordered Griffiths phase the singularity 
exponent is different from that in the disordered Grif- 
fiths phase and there is a logarithmic correction term. 
The temperature dependence of the susceptibility, which 
follows from the relation ~ fix, is given by: 



X (T) ~ T- 1+2|A| lnT, A<0 



(4.10) 



(4.8) 



We can conclude this Section that all the singularities of 
different physical quantities, both in the (strongly) or- 
dered and disordered Griffiths phases can be expressed 
by the non-linear quantum control parameter, A. 



V. DISCUSSION 

In this paper the strong disorder RG method is applied 
in the strongly disordered and strongly ordered Griffiths 
phases of the random transverse-field Ising spin chain. 
With this calculation we have demonstrated that the RG 
method leads to asymptotically exact results in the en- 
tire Griffiths region. The key concept of our solution 
is the introduction of a non-linear quantum control pa- 
rameter, A, which stays invariant under the RG trans- 
formation, even if the renormalization is approximative 
in the starting decimation steps. A, which is a mea- 
sure of the asymmetry in the renormalization between 
the couplings and the transverse fields, is simply related 
to the dynamical exponent, z, and all the singularities 
of the different physical quantities in the Griffiths phases 
can be expressed with it. In this way we have presented 
an example for an RG transformation, where the rele- 
vant non-linear scaling fieldti outside the critical fixed 
point is exactly constructed and the off-critical singular- 
ities are analytically calculated. The line of fixed points 
controlling the singular behavior in the Griffiths phases 
are found strongly attractive: for any weak initial disor- 
der, having the same asymmetry parameter, A, the sys- 
tem scales into the same fixed point. This is a remarkable 
weak-universality pijoperty of the system. We note that 
previous numericalaiij and analyticalli3 results about the 
RTIC are in accordance with our RG findings. 

At this point one may ask the question, how far these 
results are general and could apply for other random 
quantum spin systems. The above scenario is certainly 
valid for those problems, which can be mapped to the 
RTIC, so that the RG equations can be transformed 
into an equivalent form. Free fermionic spin- 1/2 mod- 
els, such as the random XX model with dimerization 
or the randon | L_0 A model with an X/Y anisotropy are 
such examplesoO. Also the one-dimensional Sinai- walk 
problem, i.e. pandom walk in a random environment with 
a global biast3, can be mapped to the RTICt2l, thus a 
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renormalization analysis^ leads to asymptotically exact 
results for this problem, too. 

Other, more general random quantum spin systems 
could belong into two main classes: i) systems having an 
IRFP and a line of semi-critical fixed points which are 
strongly attractive, and ii) those models where a cross- 
over phenomena takes place when the strength of disor- 
der is increased. Models belonging to the first class ace, 
among others, the random spin-1/2 Heisenberg chainEj 
and the random q > 2 states quantum Potts chainstil. 
For these systems we expect that a strong disorder RG- 
calculation leads to asymptotically exact results, too, and 
the physical picture obtained in the analytical treatment 
of the RTIC stays qualitatively correct. Indeed scal- 
ing arguments and numerical calculations about specific 
models are in favor of our conjecturd!3. 

Into the second class of models belong, among oth- 
ers, the random quantum clock and quantum Ashkin- 
Teller chainsE3 and probably several higher dimensional 
systems-j^two-dimensional random transverse-field Ising 
modelE3o, etc.). For these problems the disorder should 
exceed a limiting strength, when the IRFP and the line of 
semi-critical fixed points become attractive. Above this 
limiting disorder value the strong disorder RG could be 
asymptotically exact. To verify this possibility, however, 
one should perform detailed numerical investigations. 
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tt 2 (X,CI) = p 2 (A,Q)-c 2 (A) , 
and we can write a simple differential equation: 



dpi 

dr 



pi 



(A2) 



(A3) 



Since Eq. (|A3|) is e quiva lent to Eq.( 2.13| ), we have for its 
solution from Eq.( [2.16 ) with the substitution A — > c 



PI 



P° c - 



c 2 th [cm(fi A>)] 
•pfth [c\n(Cl Q /Cl)} 



(A4) 



where now p® = p(X, SIq) at a reference point, CI = Clg. 

Close to the line of fixed points, as CIq/CI — ► we should 
have c(A) — ► 0, in order to have a finite scaling combi- 
natio n in Eq . ( |A4|) . In this small A limit we obtain from 
Eqs.(|J) and that c 2 (A) = a 2 A + 0(A 2 ) where the 
value of the constant a is connected to the average lengths 
at Cl = Cl , = 7(fto). as 7 = a 2 /2(p° l ) 2 . As Cl -> the 
average length of a cluster or bond is divergent as: 



I = Iapg (ln^ 



(A5) 



which is the same as the typical distance between re- 
maining spins, Lq, as given in Eq.(2.27) if we make the 
identification, a 2 — pq. To obtain the joint distribution 
of the fields (couplings) and lengths we use the fact that 
at the fixed point of the RG transformation the appropri- 
ate scaling variable in Eq.( A.4) is aA 1 / 2 Ln(f2 /O) = O(l), 
therefore: 



Pi(X,Cl) = aA 1/2 cth 
7r z (A,ft) = aA 1/2 sh~ 1 



(A6) 



APPENDIX: DISTRIBUTION FUNCTION OF 
LENGTHS 

Here we present the general solution for the joint distri- 
bution functions, Pi(h, I, SI) and Ri(J, I, Cl), in the vicinity 
of the line of fixed points, CI/CIq — > 0. This amo unts to 
solve the set of differential equations in Eqs.( |3.9[ ) involv- 
ing the functions pi(X, Cl),m(X, Cl),ri(X, CI) and pi(X, CI), 
which appear in the Laplace-transforms, Pi(h,X,Cl) and 



Ri(h,X,Cl), in Eqs.Q) and Q . 

We start with the solution at the critical point where 
the fixed point distribution of the couplings and the 
transverse fields are identical, so that pi(X, CI) = n(A, CI) 



and 7T;(A, CI) = 
tial equations: 

dp; 

dr 



pi(X, CI). Here we have just two diffcrcn- 



-7T 



din 7T; 

dr 



-pm 



(Al) 



From Eq.(|Ai) follows that d(pf - 7r 2 )/dr = 0, therefore 
pf and 7r 2 differ by an ^-independent term: 



Consequently at the fixed point Pi(h, I, CI) and Ri(J, I, CI) 
can be obtain ed by the inv erse Laplace-transform of 
Eqs.(|J) and (|J) with Eqs.©. 

In the Griffiths region, i.e. outside the critical point 
one should consider the system of four coupled differen- 



tial equations in Eqs.(3.9), where one can construct two 



il-independent combinations of the variables: 



Pl (X,Cl)-r l (X,Cl) = 2A(X) 



and 



pi(a, ci)n(x, si) - MX, ci) Pl (x, ci) = d(x) 2 



(A7) 



(A8) 



Thus there are two variables left: yi = (pi + n)/2 and 
ui = pi — hi, which satisfy the differential equations: 



dyi 

dr 



y 2 = d 2 , d 2 = A(X) 2 + D(X) 2 , (A9) 



and 
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(A10) 



^+ U ?_( W + |A(A)|)«, = 



The solution of Eq. ( |A9| ) is analogous to that of Eq. ( |A3| 
and immediately given by: 



yi 



yfd + d 2 th [d(r - r )] 
d+yfth[d(T-T )] ' 



(All) 



where = yi (A, Op) at t he reference energy, fi = £Iq. 

To integrate Eq.flAlcj) we first notice that it is a 
Bernoulli-type differential equation and its solution can 
be expressed as: 



+ c , 



(A12) 



with 



E(T) = exp (- J [ m (T') + |A(A)|]dT') . (A13) 



Using the solution for yi(T) in Eq.(All) we can perform 
the integration for E(T) as 

E(T) = e - r l A l {d ch[d (r - r )] + y? sh[d (T - ro)]}" 1 , 

(A14) 



and putting it into Eq.(Al2) one can integrate once more 
giving: 



f + (T)-f_(T) 



f+(T)(d + |A|)-i + f-(T)(d- |A|)-i + Ce-r|A| ' 

(A15) 



where /±(r) = (yo ± d) exp[±(r — To)d] and the value of 
the constant, C, follows from the boundary condition at 

r = r . 

Having the solution at hand first, we check that at 
the critical point, where A(A) = and thus d 2 — 
D 2 = c 2 we recover the previous solution. Indeed, the 



constant in Eq.(Al5) at the critical point i s giv en by 
C = \(y?) 2 - d 2 ] L/ ' z and than combining Eq.(|All|) with 



Eq.(AlE) in the sma ll f2 and A scaling limit we recover 
the result in Eq.(A6). 

I n th e Griffiths phases, |A(A)| > 0, keeping in mind 
Eq.(3.8) we have in the small A limit: 



d(A)-|A(A)| 



d(X) 2 - A(A) 2 
d(A) + |A(A)| 



so th at th e ap propriate scaling combination in Eqs. 
dAlTl ) and ( |A15| ) as -> and A -> is A(f7 /rj) 2|A| = 



O(l). The constant, A, in Eq.(A16) is related to the av- 
erage cluster-size at Q = flo, l s (£lo), so that finally we 
obtain along the line of semi-critical points: 
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(y° + \A\) 2 (n /n) 2 ^\/(2\A\)l s (n ) 
(yf + |A|)2(f> /W |A| A(l/4A)27 s (r! ) + 1 



(A17) 



Now the joint distribution of the fields (couplings) and 
lengths can be obt ain ed through inve rse L apla ce trans- 
formation of Eqs. (|3.3| ) using Eqs. ( 3.6) and (3.7) and the 
solutions in Eqs. ( |A7|) , (As|), ( All ) and (|A17| ). 



AA|A(A = 0)| + O(A 2 ) , 
(A16) 
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